The hyperbolic Radon transform is a commonly used tool in seismic processing,for instance in seismic velocity analysis, data interpolation and for multipleremoval. A direct implementation by summation of traces with different moveoutsis computationally expensive for large data sets. In this paper we present anew method for fast computation of the hyperbolic Radon transforms. It is basedon using a log-polar sampling with which the main computational parts reduce tocomputing convolutions. This allows for fast implementations by means of FFT.In addition to the FFT operations, interpolation procedures are required forswitching between coordinates in the time-offset; Radon; and log-polar domains.Graphical Processor Units (GPUs) are suitable to use as a computationalplatform for this purpose, due to the hardware supported interpolation routinesas well as optimized routines for FFT. Performance tests show large speed-upsof the proposed algorithm. Hence, it is suitable to use in iterative methods,and we provide examples for data interpolation and multiple removal using thisapproach.
展开▼
